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TITLE OF THE INVENTION 

METHOD FOR THE FUEL-OPTIMIZED SELECTION OF A THRUSTER 
CONFIGURATION 

CROSS-REFERENCE TO RELATED APPLICATIONS 

5 This application claims priority to German Patent Application No. 102 36 570.9, 
which was filed August 8, 2002. 

FIELD OF THE INVENTION 

The present invention relates to a method for the fuel-optimized selection of a 
configuration of thrusters on a spacecraft. 

10 BACKGROUND OF THE INVENTION 

Such a method of fuel optimized selection is known for example from U.S., Patent 
No. 6,347,262 Bl for the case of a spin-stabilized spacecraft. A configuration of 
thrusters on a spacecraft, as considered by the present invention, serves in 
particular the attitude and position correction of the spacecraft. Such an 
1 5 attitude and position correction via the thrusters is known for example from EP 
0 750 239 A2. 

From EP 0 977 687 we know of a special method for the low-fuel control of an 
arrangement of thrusters on a spacecraft, wherein for the purpose of finding a 
low-fuel solution for the control of a convex linear optimization problem is 
20 resolved through 
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an initialization phase for finding a first permissible solution to the 

linear optimization problem and 

a subsequent iteration phase, in which, proceeding on the permissible 
solution to the linear optimization problem, an iterative optimization of 
an effectiveness criterion takes place. 

In this method a dual simplex algorithm is applied, which is supposed to find an 
optimal solution to the problem through a largely unfocused search method, 
wherein however it is possible with this method that there is no solution for the 
currently existing force-momental vector and the present thruster arrangement. 

From N. Karmakar: A new polynominal time algorithm for linear programming, 
Combinatorica 4 (4), 1984, p. 373-395 we know of a basic method for solving 
linear optimization problems of a general form. 

SUMMARY OF THE INVENTION 

It is an the object of the present invention to offer a method for the fuel- 
optimized selection of a configuration of thrusters on a spacecraft, which permits 
a focused search of a solution to the linear optimization problem that is 
permissible in any case. This object is achieved through the features of patent 
claim 1 described herein . 

The present invention relates to a method for the fuel-optimized selection of a 
configuration of thrusters on a spacecraft, wherein for the purpose of finding a 
fuel-optimized solution for the selection process a linear optimization problem, 
particularly a convex linear optimization problem, is resolved, through 
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- an initialization phase for finding a first permissible solution to the 
linear optimization problem and 

a subsequent iteration phase, in which proceeding on the permissible 
solution to the linear optimization problem an iterative optimization of 
an effectiveness criterion occurs. 

Pursuant to the invention it is provided that 

- in each iteration step a scaled iteration gradient is formed and 

- the iteration gradient is multiplied with a limiting factor for a 
maximum iteration threshold, which is formed while taking at least 
one boundary value condition for a permissible solution into account. 

By applying a scaled iteration gradient, instead of a mere search method as in 
the state of the art, a focused locating process for the optimal solution can take 
place. When forming the scaled iteration gradient, again especially at least one 
boundary value condition for a permissible solution ean may be included. 
Application of a scaled iteration gradient also largely excludes the circumstance 
that when only a suboptimal solution to the linear optimization problem is found. 
The fact that linear problems can involve , in particular, particularly a so-called 
convex problem is basically known, see,, for example,, the chapter "Linear 
Programming" at the following internet link of the European Business School of 
the Schloss Reichartshausen University: 

http ://w w w . ebs . de/Lehr stuehle /W irtschaftsinf or matik/NE W/Cour ses/Semester 2/ 
Math2/. 
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By taking the boundary value conditions into account within the framework of 

the limiting factor^ it is achiovod that as the rospoctivoly next iteration solution 

again a solution is may be determined. , which This solution is within a 

permissible range of values, because the limiting factor allows the iteration step 

5 width to be adjusted accordingly so that a boundary value condition is not 

violated. With a possible consideration of the boundary value conditions within 

the framework of forming the scaled iteration gradient, in particular the 

gradient direction ean may be selected substantially such that onco again a 

solution that is within a permissible range of values is determined as the 

1 0 respectively next iteration solution. 

Furthermore, for the present linear optimization problems it is known that the 
optimal permissible solution, which corresponds to a coordinate point in a 
multidimensional space of all permissible solutions that is limited by boundary 
conditions, is located on the boundary of said limited space. In this manner, the 
1 5 scaled iteration gradient and the limiting factor are preferably adjusted such 
that an iterative approximation of an optimal point on the delimitation of the 
multidimensional space of the permissible solutions occurs. 

It ean may now be provided in particular that an upper bound for a permissible 
solution is defined as a boundary value condition. 

20 Moreover it eafi may be advantageously provided that the iteration gradient is 
determined with the help of a Gauss elimination, which represents a very fast 
method. 
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It ean may also be provided in particular that in every iteration step a scaling of 

the iteration gradient takes place such that a gradient component becomes 

smaller the closer the corresponding component of the result of the previous 

iteration step comes to a boundary value condition. In this way A a new scaling 

5 operation of the iteration gradient is performed in every iteration step, wherein 

certain components of the gradient disappear when the corresponding 

components of the previous iteration solution come very close to a boundary 

value condition, irer for example,, they become smaller than a first pre-defined 

distance. Said first distance can also be selected to be infinitesimally small. 

1 0 Furthermore it ean may be provided that the iteration phase is terminated as 
soon as the result of an iteration step exceeds at least one boundary value 
condition and that the result of the previous iteration step is determined as an 
optimal solution of the effectiveness criterion. Thus the iteration is terminated if 
the algorithm leaves the range of permissible solutions, and the last permissible 

1 5 solution is determined as the optimal solution. In this way it is guaranteed in a 
simple manner that in any case the solution that is determined as the final 
result of the method is as optimal as possible and is simultaneously permissible. 

The iteration phase, however, ean may also be terminated as soon as the 
iteration method converges against a permissible solution and the result of a 
20 certain iteration step differs from the result of a previous iteration step by less 
than a second predefined distance, wherein the result of the last iteration step is 
determined as an optimal solution of the effectiveness criterion. 

DETAILED DESCRIPTION OF THE INVENTION 
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The following oxplaino one specific A preferred embodiment of the present 
invention is presented herein . 

A method for the fuel-optimized selection of an arrangement of thrusters on a 
spacecraft is considered, which is used for attitude and position control of the 
5 spacecraft. In order to generate forces and moments that are applied on a 
spacecraft, for examplei in order to be able to govern translation and rotation 
simultaneously during a docking phase or any other attitude and position 
control, n > 7 thrusters are required. The appropriate control signals must then 
meet the requirements 

10 - of being positive 

- of being smaller than a maximum value (in general equal 1). 

Furthermore, with more than 7 thrusters additionally an effectiveness criterion, 
which corresponds in general to fuel consumption, eon may be optimized. 

The mathematical formulation thus leads to the following linear optimization 
15 problem (LOP): 

Searching: thruster control a, for which the following applies: 



(a) 



0<a<f 



(b) 



T c a = r 



dim T c = 6 x n, 



(1) 



so that the following effectiveness criterion is met: 



20 



Sgi ai = g T a -> min 
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wherein 

a: control vector 

f: vector that contains the maximum values of ao 
g: vector of the weight factors 



T c = | Jli^iln | : weight and momental matrix 



r: vector that contains the required forces and moments. 

To apply all LOP solution methods, initially one permissible solution must be 
found in an initialization phase, i.e.,, a vector a z , which fulfills (la) and (lb). 
With the so-called singular value decomposition (abbreviated SVD) of T c 

10 T c = VEUi T 

with 

y. ^ orthogonal matrix 

U = (Ui, U 2 ): 
-> T c U 2 = 0 

15 E = diag al > 0: diagonal matrix of the singular values 

all solutions of (lb) can be described as 
ao = UiE-iVTr + U 2 c 
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(2) =: Ui s + v 

(2) reveals the following: 

(i) to realize random forces and moments all (al) must be greater than zero, 
i.e.^ T c must have full rank 

(ii) the first addend is completely determined by r and represents the solution 
with minimal norm of (lb) 

(hi) the second addend with the vector c that still has to be determined serves 
to fulfill the boundary condition (la) and to minimize (lc) 

From the fact that with the thruster set it must be possible to realize both 
positive and negative r, it results from (2) that ci and C2 must exist so that 

(a) ai: = Ui si + U 2 ci > 0 

(b) a 2 : = Ui (-si) + U 2 c 2 > 0 

(3) 

(c) -» U 2 (ci + c 2 ) =: U 2 c p > 0 

i.e.* the existence of a random number of vectors c p with 

(4) v p = U 2 c p > 0 

is guaranteed. After selecting a certain (a priori fixed) v p , ao e&& may be 
made positive pursuant to 

(a) ao = Ui s + v p ki > 0, 

8 
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(5) 



(b) ki = max ( Mi±£ t 8 „ 0.04 > 0 

v/>,. 



wherein e for numerical reasons was introduced for application of the following 
optimization steps. 

5 For large right sides r it is possible that (la, b) has no solution a < f , therefore 
the problem that is expanded by x s is considered 

(a) T c a = r (1 - x s ), (b) 0<a<f 

0 < x s < 1 



(6) (c) -> (To, r) 



\ x s 



= r 



10 (d) g T a + g s x s -> min 

(e) g s > g T f, so that x s becomes zero if r is not too large 

This now also allows the upper bound to be adhered to and allows the required 
permissible starting value for a to be calculated for 

(a) a z = ao (1 - x sz ) = Ui s (1 - x sz ) + v p ki (1 - x S2 ) 
15 ^ w, w > 0 

(7) (b) Xsz = 

0, w < 0 
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(c) 



w = max < 1 - 



► 



s here represents a first, for example^ infinitesimal, distance. 

All subsequent considerations relate to the expanded system (6), wherein 
however the original description pursuant to (1) is maintained for reasons of 
simplicity. 

To resolve the LOP a second procedural step now follows, namely an 
optimization of the efficiency criterion (1c) and/or (6d), which is performed 
iteratively as follows: 



In this Vgi represents the iteration gradient, which is scaled in every 
iteration step, i.e^ in each iteration step the gradient direction is newly 
determined. Moreover, k represents a limiting factor for the iteration step 
width, which is determined as follows: 



ai+i = ai - Di U 2 < i > U 2 ® T Di g • k 



(a) 



vgi • k 



(b) k = min (k u , ko) (1 - e) 



(8) 





~fj 




(c) 



mm 



V 
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This selection of k while taking the boundary value condition 0 < a < f into 



account ensures that ai+i remains permissible. 



(d) 



D t = diag a 




,y = l,..,n + l 



(e) U 2 ®: zero space vectors of T c Di = Vd) E® UiO) T 



5 



ai = a z , calculated pursuant to (7). 



The essential idea in (8) is the constant scaling of the problem with Di and the 
subsequent continuation into the thus modified, on U2 (i) projected, negative 
gradient direction, wherein, due to the familiar structure of the problem as a 
convex linear optimization problem with an optimal solution on the boundary, it 
10 is guaranteed that the efficiency criterion is reduced in every iteration step. The 
iteration is preferably interrupted when the amount of Vgi • k drops below a 
specified threshold as a second distance, i.e. ± ai hardly changes any more. 

A particular expansion of the present method as compared with the method ftt in 
Karmakar consists in taking an additional boundary value problem into account 

1 5 with every iteration step, here the inclusion of the upper bound f (upper bound 
problem) by adding the second factor in Di and taking the upper bound f into 
consideration within the framework of the term ko in the calculation of k 
pursuant to (8b). Up to now, in the Karmakar methods usually complex 
expansions of the linear optimization problem with slack variables were offered 

20 with the disadvantage that thon the dimension of the problem that needs to be 
solved is increased considerably. Here the present method represents an 
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essential simplification. Additionally it eliminates the very complex 

determination of a permissible solution in the initialization phase, as practiced 

in Karmakar, through the suggested initialization phase, which is better adapted 

to the present problem. 

5 Another advantageous procedural step of the method described here is thus in 
the constant calculation of Vgi, which preferably occurs not through the SVD of 
T c Di, but due to (here we use the following simplified depiction: T c for T c Di, D for 
DO 

(a) M = T c Tc T = V E 2 Vr 
10 -> 

(b) UiUi T = Tc T M-i T c 

v g = D U 2 U 2 T Dg = D(I - UiUi T ) Dg 

b 

(9) = D (I - Tc T M-i T c ) b 

(c) 

15 =D(b-T c T x ) 
with 

(d) Mx = Tcb 

via the solution from (9d) for x by means of a Gauss elimination. This method is 
clearly faster than an SVD. For a 0 as well, Ui is not determined, but Ui s is 
20 calculated directly pursuant to 
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(e) 



Uis = Ui S-iVT r = Tc T M- 1 r 



Thus the method consists in the execution of the calculation steps 



- equation (5) with (9e) 



} 



- equation (7) 



initialization 



equation (8) with (9c, d) 



} 



iteration 



Finally additional advantages of the suggested method compared with a simplex 
method pursuant to the state of the art are summarized in the following: 

- a calculation of a permissible solution can occur with minimal effort 
(even a z could already be used as a thruster selection, although it 
would not be fuel-optimized), 

the subsequent optimization can be limited to a few steps if necessary 
when* e.g.* the power of the on-board computer is limited, 

experience shows that the optimum on average is achieved with 
considerably fewer arithmetic operations and generally requires less 
memory than with a simplex method, 

- in the case of a thruster failure, T c must be reduced only by the 
appropriate column, and an appropriate vector v p must be stored, i.e.* 
in this case only minimal additional memory is required, in contrast to 
the state of the art, 



Attorney Docket: 010408.52553US 
MARKED -UP VERSION OF SPECIFICATION 

by including the upper bound problem in practice problems are avoided 

especially with limited thruster adjusting capacity. 
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